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1.  INTRODUCTION 

Cross-classified  data  having  ordered  categories  arise  in  many  investiga¬ 
tions  conducted  by  medical ^  physical,  natural  and  social  scientists.  Statistical 
methods  have  been  developed  and  continue  to  evolve  to  analyze  such  data.  Many 
of  these  methods  are  tailored  to  answer  specific  questions  and  Issues  raised. 

For  review*  of  the  literature  in  this  area,  the  reader  is  referred  to  Agresti  [1] 
and  Goodman  [6]. 

We  begin  with  a' general  description  of  a  problem  tackled  in  this  paper. 

Let  8  be  the  Borel  a  -  field  on  the  real  line,  R  and8*B  the  product  a  -  field 
on  RxR.  Let  y  be  a  probability  measure  on  8^8  and  yj^  and  y2  the  corresponding 
marginal  probability  measures  on  8  ,  i.e.,  y^B)  =  y(BxR)  and  y2(B)  *  y(RxB) 
for  every  B  in  B,  Following  Lehmann  [14],  y  is  said  to  be  a  positive  quadrant 


dependent  i  f 


u{[c  ,«)  x  [d  _>  p  |  [c  ,-)>  \x  \  [d  ,»)/ 

for  every  c,d  in  R.  In  the  jargon  of  random  vairables,  the 

above  notion  can  be  rephrased  as  follows.  Let  X  and  Y  be  two 

random  variables  with  some  joint  probability  distribution 

function  F.  X  and  Y  are  said  to  be  positive  quadrant  dependent 

P{X  >  c,  Y  >  d}  >  Pjx  >  c}  P{Y  >  d> 

for  all  c,d  in  R.  For  various  properties  of  positive  quadrant 

dependence,  see  Lehmannjj^)  or  Eaton  [3].  In  this  paper,  we 

look  at  the  notion  of  positive  quadrant  dependence  from  a  global 

point  of  view.  Let  M  denote  the  set  of  all  positive 

PQD 

quadrant  dependent  probability  measures  u  on  IB  xfl  .  It  is 

natural  to  think  along  the  following  lines.  If  M  is  a  convex 

PQD 

set  and  compact  in  some  decent  topology,  then  the  set  of  extreme 

points  of  M  will  be  non-empty.  See  Phelps  [Jcjj .  Moreover, 

PQD 

every  member  of  M  can  be  expressed  as  a  mixture  (in  some 

PQD 

sense)  ot  extreme  points  of  M  .  There  are  certain  properties 

PQD 

of  d i s t r i bu t i ons  which  are  preserved  under  mixtures.  Under  these 

circumstances,  it  suffices  to  examine  extensively  the  extreme 

points  so  as  to  make  comments  on  the  members  of  M  .  But  this 

PQD 

line  of  reasoning  fails  since  M  is  not  a  convex  set  as  the 

PQD 

the  following  example  demonstrates. 


Let  y  be  a  probability  measure  on  (&  x<B  with 
support  contained  in  {(1,1),  (1,2),  (2,1),  (2,2)}.  Such  a 


probability  measure  can  be  written  as 
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if  and  only  if 
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We  can  identify  some  natural  subsets  of  M  as 

PQD  ra 

convex  sets.  Let  x  and  v  be  two  probability  measures  on  UJ  . 

Let  M  (  x  ,  v )  be  the  collection  of  all  probability  measures  y 
PQD 

in  M  such  that  y  =X  and  y  =  v ,  i  .  e  .  , 

PQD  1  2 

M  (x,v)  =  { y  e  M  ;  y  =  X  and  y  =  vl. 

PQD  PQD  1  2 

In  Section  2 ,  we  show  that  M  (  x  ,  v )  is  a  compact  convex  set 

PQD 

in  the  weak  topology  on  the  space  of  all  probability  measures 

M  on  (B  x®  .Using  this  result,  one  obtains  a  decomposition  of 

M  as 
PQD 

M  =  \J  U  M  (  X  ,  v )  , 
pqo  X  v  PQD 

where  the  union  is  taken  over  all  probability  measures  x  ,  v 
on  &  . 

In  Section  3,  we  concentrate  on  the  case  when  both 

>  and  v  have  finite  support.  We  describe  a  method  of 

enumerating  all  extreme  points  of  M  (  x ,  v)  with  the  help  of 

PQD 

some  examples.  In  Section  4,  using  the  structure  of  M  (x,v), 

PQD 

we  compare  the  performance  of  some  tests  for  testing  independence 
against  strict  positive  quadrant  dependence. 


2 .  Main  Resu  1  ts 

In  this  section,  we  show  that  for  any  two  probability 

measures  x  and  v  on  &  ,  M  ( x ,  v )  is  compact  and  convex. 

PQD 


•  «*  *•  -■  V  .* 
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We  need  the  following  definitions  and  results  in  this  connection. 


Let  (X,d)  be  a  Polish  space,  i.e.,  a  complete 
separable  metric  space.  Let  QJ  be  the  Bore!  a-field  on 


X 

and  the  space  of  all  probability  measures  on 

equipped  with  weak  topology. 


X 

i  s 


Definition  1 .  A  subset  S  of  M  is  said  to  be  uni formly 

A 

tight  if  for  every  e  >0  there  exists  a  compact  subset  C  of 
X  such  that 

y  ( C )  >  1  -  e  for  every  y  in  S. 


The  following  is  known  as  Prohorov's  theorem. 

Proposition  2  A  subset  S  of  M  is  relatively  compact  if 

X 

and  only  if  S  is  uniformly  tight.  S  is  compact  if  and  only  if 
S  is  closed  and  uniformly  tight. 

Proof.  See  Billingsley  [.  ,  Theorems  6.1  and  6.2,  p.37]. 


Theorem  3  Let  M(x,v)  be  the  collection  of  all  probability 

measures  y  on  (B  x®  such  that  y  =  x  and  u  =  v  . 

1  2 

Then  M(x,v)  is  compact. 

Proof .  It  is  obvious  that  M(x,v)  is  a  closed  subset  of  M  , 
the  space  of  all  probability  measures  on  ®  x(B  .  We  show  that 
M ( x , v )  is  uniformly  tight.  Let  e  >  0.  There  exist  compact 
subsets  C  and  C  of  R  such  that  x(Cc)  <  e/2  and 


1  2  l 

v(Cc)  <  e/2.  C  xC  is  a  compact  subset  of  RxR.  Let  y  e  M(x,v). 

2  12 

Then 


-6- 

y  C(C1xC2)C^  <  y (C^xR  U  RxC  ^  ) 

<  y(cJxR)  +  y ( RxC ^ ) 

=  UjtcJ)  +  y2(Cp 
=  x(cj)  +  v(C2) 

<  e 

This  completes  the  proof  in  view  of  Proposition  2  . 

The  following  result  is  the  main  result  of  this  section 
Theorem  4  For  any  given  probability  measures  X  and  v  on  Ti, 
Mrq0(a,v)  is  compact  and  convex. 

Proof .  MpQD(x,v)  is  a  closed  subset  of  M(x,v)  follows  from 
the  following  observation.  Let  yn,  n>l  be  a  sequence  in 
MpQo(*>v)  converging  weakly  to  a  y  in  M(x,v).  Then  for  any 
c,d  in  R,  (See  Billingsley  |'2»  P-  11  I  )  » 

u  {  !"c  » ®  )  x  | 'd  ,  ® ) }  _>  lim  sup  yn{|c,°Ox|d,®)} 

n-*-°° 

_>  X  {  |  c  ,  ® ) }  v  {  f  d  , » ) }  . 

Hence  y  e  MpQD(x,v).  This  implies  that  MpQD(x,v)  is  compact. 
We  now  show  that  M  p  ^  D  (  X  ,  v  )  is  convex.  Let  y,n  e  MpQD(\,v) 
and  0  <  a  <  1.  Then  for  any  c,d  in  R, 


_■  -a 


(au  +  ( 1  -  a  )  n  )  {  [c  , «  )  x  [d  , » ) }  -  x{[c,»)}v{fd,^)} 

=  a  p  {  [c  ,  ®  )  X  [d  ,  00 )  }  +  ( 1  -  a  )  n  (  [c  ,  °° )  x  [d  ,  ®  )  }  - 

x{  [c  , « ) } v {  [d  , o° ) } 

±  ctX{[c,oo)}  v{Qj,°°)}  +  ( 1  -  a)X{  [c,oo)}v{  [d  ,°o)} 

X  {  [c  ,00  )  }v  {  [*d  ,0=')  }  =  0  . 


Consequently,  ay  +  ( 1  -  a  )  n  e  MDnn(x,v).  This  completes  the 


3.  Extreme  Points 


In  this  section,  we  assume  that  the  support  of  a  is  {l,2,*»*,m} 

and  that  of  v  is  {1,2, Let  p..  =  u({i}),  i  =  l,2,***,m  and  q.  = 

J 

j  =  l,2,*»*,n.  In  this  case,  we  use  the  suggestive  notation 

MPQD^p1  ,p2  *  ,pm’  ql,q2’"*’qn^  for  Mpqd^a*v^  Any  v  in  mPq0(^)  can 

be  written  in  the  following  form 


pll 

P12  *•* 

pln 

rH 

CL 

P21 

P22 

p2n  ! 

P2 

• 

•  •  • 

•  •  • 

•  •  • 

• 

• 

pml 

Pm2  *** 

pmn 

Pm 

ql 

<^2 

qn 

1 

where  p^.  =u((i,j)},  i  =  l,2,«*»,m;  j  =  l,2,»**,n.  In  other  words, 

MPQD^p1’p2’”  *  ,pm’  ql,q2’’”,qn^  is  the  co11ection  °f  a11  matrices  (p^.) 
of  order  mxn  such  that  each  p„  j>  0,  row  sums  Pi  »P2»*  *  * » Pm »  column  sums 
ql,q2’**"’qn  and  the  Joint  distribution  is  positive  quadrant  dependent.  The 
compact  convex  set  ^pqq(  Pi  »Pj  *  ’  *  *  >Pm»  ql,q2’”‘,qn^  ^as  a  finite  number  of 
extreme  points.  We  now  describe  a  method  of  enumerating  the  extreme  points 
in  some  special  cases  for  illustration  from  which  the  general  technique  can 
easily  be  perceived.  As  we  shall  see  shortly  in  Section  4,  the  knowledge 
of  extreme  points  has  considerable  bearing  on  the  power  of  tests  of  indepen¬ 
dence  aoainst  strict  positive  quadrant  dependence. 
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Example  1.  m  =  2  and  n  =  2. 

Let  Pj  *P£ »  ql’q2  be  sPec1^ied  such  that  each  °f  Pl*P2’ql,q2  1S 
positive  and  Pj  +  P2  =  1  =  Pj  +  P2*  It:  can  156  easi^y  verified  that  if  a 
matrix 


"i’ll 

P12 

.P21 

P22 

row  sums 

P 

t°  MpnD(Pi>P2;  qi*q2)’  then 


p2q2  -  p22  -  p2  A  q2 ' 


Conversely,  if  the  number  p22  satisfies  the  above  inequality,  then  the  matrix 


prVp22 

d2-P 

P2-P22 

P22 

belongs  to  mpqq ( Pj  » P2 ;  ^1  »^2 ^  *  T*iere  are  two  extreme  points  of 
MPQd( P1  ,P2 ;  ql»q2^*  These  are  given  by 


plql 

plq2 

p2ql 

p2q2 

plql 

plq2 

P2qi 

p2q2 

ql  q2"P2 
0  P, 


if  p2  A  q2  =  p2. 


P2'q2  q2 


if  p2  A  q2  =  q2 , 


Every  member  of  P2  *  p2  *  qi,q2^  is  a  convex  combination  of  these  two 

extreme  points. 


v  \  *.  *.  . 


'  .  ‘  k  *  '  k  .  *  ■>  V  ?.  -  * 
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Example  2  m  *  2  and  n  =  3. 

Let  Pi *P2 ,q2 ,q3  ^ve  P0S1tlve  numbers  satisfying  +  p2  = 

1  =  qx  +  q2  +  q3 .  If  the  matrix  (p. ^ )  belongs  to  mPqD(Pi»P2;  then 

p2q3  -  p23  -  p2  A  q3  and 

(p2q2  +  p2q3)  v  P23  1  P22  +  P23  -  p2  A  (q2  +  p23^' 

where  a  y  b  indicates  the  maximum  of  the  numbers  a  and  b.  Conversely, 
if  p22  and  p23  are  two  numbers  satisfying  the  above  inequalities,  then 
the  matrix 

ql'p2+p22+p23  q2'P22  q3~p23 

p2'p22"p23  p22  p23 

belongs  to  mpqd^Pi*p2;  ql ,q2’q3^  •  T>1e  i,nPact  of  this  observation  is  that 
the  numbers  p22  and  p23  in  the  matrix  ( p.  ^ )  determine  whether  the  matrix 
(p^.)  belongs  to  mpqq(Pi»P2;  *q2 * q3 )  or  not-  These  two  inequalities 
determine  a  simplex  in  the  p22  -  p23  plane.  As  a  simple  illustration,  take 
pl  =  p2  =  h\  q-^  =  q2  =  q3  =  1/3.  The  determining  inequalities  are 

1/6  <_  p23  £  1/3  and 

1/3  =  p23  v  1/3  <  P22  +  P23  £  h  A  (1/3  +  p23)  =  4- 
These  inequalities  determine  the  following  simplex  in  the  p22  -  p23  plane. 


There  are  four  extreme  points  of  the  set  mpqq( 1/2 ,1/2 ;  1/3, 1/3, 1/3) 
given  by 


1/6 

1/6 

1/6  ~ 

1/6 

1/3 

1/6 

1/6 

1/6  J  , 

_  !/6 

0 

1/3  1/6 


0  1/6 
1/3  1/6 


L  0  1/6  1/3  J  ,  [_  0  1/3  1/6 | 

corresponding  to  the  four  points  P^,  P^,  P^,  P^  respectively.  Every 
member  of  MpQg(l/2,l/2;  1/3, 1/3, 1/3)  is  a  convex  combination  of  these 
four  extreme  points. 


2- 


Example  3  m  =  3  and  n  =  4. 

Let  pi>p2»P3»qi»q2»q3»q4  be  seven  positive  numbers  satisfying 

P1  +  p2  +  p3  =  1  =  ql  +  q2  +  q3  +  V  If  ^pij-)  e  MPQD(p1’P2’P3;  ql *q2’q3*q4^ 
then 

O)  P3q4  i  P34  1  P3  A  q4  * 

(2)  P34  V  (P3q3+P3q4)  1  P33  +  P34  1  P3A(q3+P34)> 

(3)  P34  V  (p2q4+p3q4)  <  p24  +  p34  <  q4  A  (p2+p34), 

(4)  (p33+p34}  V  {p3q2+p3q3+p3q4)  ±  p32  +  p33  +  p34  1 

P3  A  (q2+P33+P34). 

(5)  (P33+P24+P34)  V  (p2q3+p2q4+p3q3+p3q4)  <  p^  +  P24  +  P33  +  P34 

-  ^ P2+P33+P34^  A  ^q3+p24+p34^  ’ 

and 

(6)  (P23tP24+P34+P32*P33)  *  <p2  VP2VP2  i 

P59  +  P?-3  +  P9/1  +  P'39  +  P'3'3  +  P34  1 


(p2+p32+p33+p34^  A  ^q2+p23+p24+p33+p34^ 


Conversely,  if  p34,  p33,  p24,  P32,  p^,  p22  are  six  numbers  satisfying  the 
above  inequalities,  then  these  six  numbers  determine  uniquely  a  member  of 
MPQD^p1,p2’p3’  ql ,q2,q3 ,q4^ '  We  exP^aln  b°w  to  enumerate  all  extreme  points 
of  the  convex  set  in  the  simple  example  Pj  =  P2  =  p3  =  ^  and  q^  =  q2  = 
q3  =  q4  =  1/4.  (The  technique  in  the  general  case  is  similar  to  this  special 
example.)  The  six  inequalities  above  now  become 


j r*y*.)F\*\^J*y*\*\'*\,*V*\YV*  '.ip  E*?W gTS*r»^»  L"*.T».^l^V-  L^U.-g IT* l-  V-»>-» ir» V-» y- >r^ ir- >r^ >  - 1 
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(1) 

1/12 

±  p34 

£  3/12, 

(2) 

p34 

V  2/12  j 

i  P33 

+  P34 

< 

4/12  A  (3/12 

+  p34> 

=  4/12, 

(3) 

p34 

V  2/12  . 

1  p24 

+  P34 

< 

3/12  A  (4/12 

*  p34> 

=  3/12, 

(4) 

(p33 

+  p34>  ’ 

V  3/12 

1  P32  + 

P33  +  P34  - 

4/12  A 

(3/12  +  p33  +  p34) 

=  4/12, 

(5) 

<p33 

+  p24  + 

P34)  ’ 

V  4/12 

< 

p23  +  p24  +  » 

J33  +  p34  - 

(4/12 

+  P33 

+  P34 ' 

1  A  (3/12  +  p.,4  + 

p34>  * 

3/12  +  p^4  +  P34» 

(6) 

(P23 

+  P32  + 

P33  + 

P24  + 

p34> 

V  6/12  £ 

P22  + 

p23  + 

p24  + 

P32 

+  p33  +  p34 

< 

(4/12 

+  P32 

+  P33 

+  P34)  A  (3/12  + 

p23  +  P 

24  +  p33  +  p34^  * 

The  first  step  in  the  determination  of  extreme  points  is  to  get 


rid  of  the  maximum  and  minimum  symbols  by  splitting  some  or  all  inequalities 
above.  For  example,  inequality  (1)  can  be  written  as  1/12  <  p34  £  2/12 
and  2/12  £  p34  £  3/12.  Inequality  (3)  can  be  written  as  p34  V  2/12  £ 


p33  +  p34  i 

3/12 

and 

3/12 

< 

P33  +  p34 

£  4/12.  The 

above  set  of 

inequalities 

are  equivalent  to 

the 

following 

four  sets  of 

inequalities 

I  (1) 

1/12 

< 

p34  ^ 

2/12 

(2) 

2/12 

< 

p33  + 

p34 

£  3/12 

(3) 

2/12 

< 

p24  + 

p34 

£  3/12 

(4) 

3/12 

< 

P32  + 

P33 

+  p34  1 

4/12 

(5) 

Same 

as 

above 

(6) 

Same 

as 

above 
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The  maximum  and  minimum  symbols  in  inequalities  (5)  and  (6)  stay 
put  in  spite  of  splitting  the  inequalities  (1)  and  (2).  In  order  to  neutralize 
these  symbols,  we  introduce  the  following  auxiliary  inequalities. 

(2,3)  p^^  +  P24  P34  £  4/12 


4/12  <  p33  +  p24  +  p34 

(2, 3, 4, 5)  p23  +  P 32  P^3  ^  P24  *  P34  1  6/12 


6/12  <  p23  +  p32  +  p33  +  p24  +  p34 

(4,5)  4/12  +  p32  +  p33  +  p34  <  3/12  +  p23  +  p24  +  p33  +  p34 


3/12  +  p23  +  p24  +  p33  +  p34  -  4/12  +  p32  +  P33  +  p34‘ 

Now,  a  choice  of  each  of  the  auxiliary  inequalities  (2,3),  (2, 3, 4, 5) 
and  (4,5)  is  appended  to  each  set  of  the  inequalities  I,  II,  III  and  IV.  This 
would  generate  32  sets  of  inequalities  equivalent  to  the  four  sets  I,  II,  III 
and  IV  of  inequalities.  To  save  space,  we  will  not  reproduce  these  32  sets 
of  inequalities.  A  sample  set  of  inequalities  is  produced  below  for  further 


discussion. 
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A  sample  set  of  inequalities  chosen  from  32  sets  of  inequl ities  descri bed  above 


(1) 

2/12 

1  p34 

<  3/12 

(2) 

3/12 

±  p33 

+  P34 

<  4/12 

(3) 

''  P34 

i  P24 

+  p34 

<  3/12 

(2,3) 

P33  + 

P24  + 

P34  i 

4/12 

(4) 

P33  + 

p34  1 

P32 

+  p33  +  p34  1  4/12 

(2, 3, 4, 5) 

6/12 

±  P23 

+  P32 

+  p33  +  p24  +  p34 

(5) 

4/12 

1  p23 

+  P24 

+  p33  +  p34  -  3/12 

+  P; 

24  +  p34 

(4,5) 

4/12 

+  p32  + 

p33  + 

P34  £  3/12  +  p23  + 

p24 

+  p33  +  p34 

(6) 

P23  + 

p32  + 

P33  + 

p24  +  p34  1  p22  +  p 

23  + 

P24  +  p32  + 

<  4/12  +  p32  +  p33  +  p34. 


The  above  set  of  inequalities  is  obtained  from  the  set  IV  of 
inequalities  by  appending  the  first  choice  of  (2,3),  the  second  choice  of 
(2, 3, 4, 5)  and  the  first  choice  of  (4,5). 

In  order  to  obtain  a  member  of  MpgD ( 1/3 ,1/3, 1/3;  1/4, 1/4, 1/4, 1/4) 
we  proceed  as  follows.  Set  the  central  expression  in  each  of  the  main  six 
inequalities  equal  to  the  quantity  either  on  the  left  or  the  right  of  the 
inequalities  and  then  solve  the  system  of  equations  thus  arize  in  the  unknowns 

p34’  p33’  p24’  P32’  p23’  P22  makin9  su-'e  that  the  constraints  imposed  by 
the  auxiliary  inequalities  are  satisfied .These  system  of  equations  are  easy 
to  solve.  The  solution  will  give  a  member  of  mpqq( 1/3.1/3.1/3;  1/4, 1/4, 1/4, 1/4) 


V  V.'-'AV.T' V.V.1 A 
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Generate  members  of  MpQ[)(l/3,l/3,l/3;  1/4, 1/4, 1/4, 1/4)  by  following 
the  above  procedure  for  each  set  of  the  32  sets  of  inequalities.  The  set  of 
extreme  points  of  the  convex  set  of  interest  is  a  subset  of  these  solutions. 
There  will  be  a  large  amount  of  duplicates  and  some  of  the  solutions  obtained 
are  already  convex  combinations  of  other  solutions.  After  considerable  amount 
of  weeding,  we  got  the  following  matrices  as  the  entire  collection  of  extreme 
points  of  the  convex  set  Mp^D( l/3»l/3,l/3;  1/4, 1/4, 1/4, 1/4) . 


1. 

1/12 

1/12 

1/12 

1/12 

2. 

2/12 

0 

1/12 

1/12 

1/12 

1/12 

1/12 

1/12 

0 

2/12 

1/12 

1/12 

_1/12 

1/12 

1/12 

1  / 1 2_ 

1/12 

1/12 

1/12 

1/12  . 

3. 

ri/l  2 

1/12 

1/12 

1/12" 

4. 

"2/12 

0 

1/12 

1/12' 

2/12 

0 

1/12 

1/12 

1/12 

1/12 

1/12 

1/12 

_  o 

2/12 

1/12 

1/12. 

_  0 

2/12 

1/12 

1/12. 

5. 

'3/12 

0 

0 

1/12' 

6. 

1/12 

2/12 

0 

1/12' 

0 

1/12 

2/12 

1/12 

1/12 

0 

2/12 

1/12 

_  0 

2/12 

1/12 

1/12_ 

J/12 

1/12 

1/12 

1/12 

7. 

'2/12 

1/12 

0 

1/12“ 

8. 

"3/12 

1/12 

0 

0 

0 

1/12 

2/12 

1/12 

0 

0 

2/12 

2/12 

1/12 

1/12 

1/12 

1/12  j 

0 

2/12 

1/12 

1/12. 

9. 

”2/12 

0 

2/12 

0 

10. 

1/12 

1/12 

2/12 

0 

1/12 

1/12 

0 

2/12 

2/12 

0 

0 

2/12 

_  0 

2/12 

1/12 

1/12_ 

_  0 

2/12 

1/12 

1/12. 

2 

2 


0 

3/12 

1/12_ 

0 

2/12 

1/12' 

2/12 

1/12 

0 

1/12 

0 

2/12 

1/12 

0 

1/12' 

1/12 

3/12 

0 

1/12 

0 

2/12 

0 

1/12 

1/12 

1/12 

2/12 

0 

2/12 

0 

2/12. 

2/12 

1/12 

0 

0 

2/12 

1/12 

1/12 

0 

2/12 

2/12 

0 

0 

0 

3/12 

1/12 

1/12 

0 

2/12 

1/12 

2/12 

0 

0 

1/12 

1/12 

II  * 

■■ 

H 


1/12 

1/12 

0 

2/12 

1/12 

1/12 

1/12 

1/12 

2/12 

0 

2/12 

0 

0 

2/12 

0 

2/12 

2/12 

1/12 

0 

1/12 

1/12 

0 

3/12 

0 

0 

2/12 

0 

2/12 

2/12 

1/12 

1/12 

0 

0 

1/12 

2/12 

1/12 

1/12 

1/12 

0 

2/12 

2/12 

0 

2/12 

0 

1/12 

1/12 

1/12 

1/12 

.  0 

2/12 

0 

2/12 

1/12 

1/12 

1/12 

1/12 

1/12 

1/12 

2/12 

0 

!  1/12  2/12  0 

0  1/12  3/12  . 


4.  Testina  independence  against  strict  positive  quadrant  dependence 


Let  X  and  Y  be  two  random  variables  with  known  marginal 
distributions  and  unknown  joint  distribution.  We  want  to  test  the  hypothesis 
that  X  and  Y  are  independent  against  the  hypothesis  that  they  are  strictly  posi 

tive  quadrant  dependent.  By  strict  positive  quadrant  dependence  we  mean  positive 
quadrant  dependence  but  not  independence.  The  data  consist  of  N  independent 
realizations  of  the  vector  (X,Y).  Let  x  be  a  test  proposed  for  testing 
the  hyposthesis  of  independence  based  on  the  given  data.  Let  a  be  the 
distribution  of  X  and  v  that  of  Y.  Let  mpqq(*>v)  be  the  collection 
of  all  bivariate  distributions  with  fixed  marginals  X  and  v  which  are 
positive  quadrant  dependent.  The  power  function  of  the  test  x  can  be 
defined  formally  as  follows. 


8  (vi)  =  Pr{  x  rejects  the  null  hypothesis  /  u) 


for  u  in  M  nn(x,v).  The  above  probability  is  computed  when  the  joint 


distribution  of  X  and  Y  is  y .  The  calculation  of  the  power  function 

of  the  test  t  whose  domain  of  definition  is  mpqq(^»v)  is  very  tedious. 

Moreover,  if  we  wish  to  compare  the  performance  of  two  tests  to  discriminate 

the  null  hypothesis  of  independence  against  the  alternative  hypothesis  of 

strict  positive  quadrant  dependence,  we  need  to  compare  their  power  functions 

This  comparison  then  becomes  doubly  more  difficult  to  achieve.  But  the 

following  theorem  asserts  that  it  suffices  to  compare  the  powers  at  the 

extreme  points  of  mPqDU.v)  only. 

12  k 

Theorem  5  Let  y,  y  ,  y  y  be  members  of  Mpg0(x,v)  such 

12  k 

that  y  =  cijy  +  o^u  +  •••  +  a^y  for  some  ,  a2’  ’’’  ’  a  k  —  ® 

k 

with  z  a.  =  1.  Then 
1=1  1 

k 

B  (u)  =  £  a .  6  (y 1 ) 

T  i=l  1  T 

The  above  result  can  be  used  as  follows.  Suppose  each  of  X  and  Y  takes 
finitely  many  values.  Then  the  convex  set  MpqD(x,v)  has  finitely  many 
extreme  points  and  every  member  of  Mpgg(x,v)  can  be  written  as  a  convex 
combination  of  these  extreme  points.  Then  the  power  of  the  test  x 
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evaluated  at  any  given  u  in  Mp^(x.v)  is  precisely  the  same  convex 
combination  of  the  powers  of  the  test  evaluated  at  each  of  the  extreme  points. 
This  result  also  points  out  that  in  order  to  compare  the  performance  of 
two  given  tests,  it  suffices  to  compare  the  powers  evaluated  at  the  extreme 
points.  As  an  illustration,  we  consider  the  case  when  X  takes  values  1  and  2, 
and  Y  takes  values  1,  2  and  3.  Let  n..  .  =  total  number  of  (X,Y)'s 
with  X  =  i  and  Y  =  j,  i  =  1,2  and  j  =  1,2,3.  The  data  can  be  arranged 
in  the  form  of  a  contingency  table  as  follows. 


11 

n12 

n13 

21 

n22 

n23 

In  this  section,  we  compare  the  performance  of  two  tests  for  testing 
the  null  hypothesis  of  independence  against  the  alternative  of  strict  positive 
quadrant  dependence  in  the  context  of  2x3  contingency  tables  above. 


Tj  i  Test  based-on  gamma  rati< 


Let  the  bivariate  distribution  of  X  and  Y  be  given  by 


11 

P12 

P13 

21 

p22 

P23_ 

The  Gamma  Ratio  (see  Goodman  and  Kruskal  T25])  of  X  and  Y  is  defined  by 

IT  - 

y  _  c _ d 


n  ,  it  j 

c  +  d 


where 


2pll(p22  +  p23}  +  2p12p13  and 


-C  -’-LL-Zd  rdi'  --lZ'lJ  nd  “  2p13(p21  +  p22)  +  2  p  1 2  p2 1 

One  can  show  that  y  =  0  if  X  and  Y  are  independent.  See 


>  0  if  X  and  Y  are 


Agresti  (1984,  p.160).  One  can  also  show  that  y 
positive  quadrant  dependent.  An  estimate  of  y  based  on  the  sample  given 
above  is  given  by 

y  _  C  -  0 
C  +  D  ’ 


where  C  =  The  total  number  of  concordant  pairs  =  n^(n22  +  n23^  +  n12n23 
and  D  =  The  total  number  of  discordant  pairs  =  ni 3 ^ n21  +  n22^  +  nl 2n21 * 

The  following  is  a  natural  test  based  on  y  for  testing  the  above  null 
hypothesis  against  the  specific  alternative  mentioned  thereby. 


Test  Tj  :  Reject  the  null  hypothesis  if  y  >  a. 


T2  :  Test  based  on  eigen  values 


Let  the  marginal  distribution  of  X  be  given  by  p,  and  p2>  and 
that  of  Y  by  q1 ,  q2  and  q^.  Let 


Pll 

P12 

p13  ” 

^TT 

^lq2 

^lq3 

Q  = 

P21 

P22 

P23 

^2ql 

*^2q3 

Let  Kj  and  <2  be  the  eigen  values  of  QQ^,  where  T  denotes 
operation  transpose  on  matrices.  We  give  some  properties  of  these  eigen  values 
below.  For  further  details,  see  Lancaster  [11]  and  [12],  O'Neill  ([17],  [18], 


Properties 

1.  One  of  the  eigen  values  is  always  equals  unity.  Let  us  use 
for  this  eigen  value. 

2.  If  X  and  Y  are  independent,  ~  0* 

3.  If  X  and  Y  are  strictly  positive  quadrant  dependent,  then  >  0. 


We  estimate  based  on  the  data  given  above  as  follows 

Let 


B 


n12  n13 

N  vPpf 

n22  n23 

N/prjrrjT 


Let  and 

as  an  estimator  of 


be  the  eigen  values  of  BBT. 
k  ^  k  2  * 


Then  we  propose  <  =  < 


Test  T^  :  Reject  the  null  hypothesis  if  <  >_  a. 

We  discuss  the  performance  of  these  two  tests  in  the  case  of  two 
specific  examples  given  below. 


Example  1  P1  =  p2  =  ^  an<^  ql  =  p2  =  ^3  = 


Example  2 


p,  =  1/4,  p,  =  3/4  and  q.  =  1/2,  q,  =  1/4,  q.  =  1/4. 
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Let  us  now  elaborate  on  some  of  the  properties  of  the  eigen 
T  T 

values  of  QQ  .  The  second  eigen  value  of  QQ  ,  <2»  catl  w°rked  out 

explicitly. 

(pllp22  ~  P12P21)  +  (pllp23  ~  p13p21)  + 

K2  =  plp2qlq2  +  plp2qlq3 

(p12p23  ~  P13P22) 
plP2q2q3 

From  this  it  follows  that  =  0  if  and  only  if  X  and  Y  are  independent 


At  this  juncture,  we  want  to  make  some  remarks  on  the  definition 


of  the  matrix  B  above.  In  order  to  develop  an  estimator  of  x^,  it  is 


natural  to  divide  each  frequency  n^  in  B  by  the  square  root  of  the 


product  of  the  corresponding  marginal  totals  n^  and  n^ .  See  O'Neill 
([17],  [18],  [19]).  If  we  had  proceeded  as  outlined  by  O'Neill,  one  of  the 


eigen  values  of  BB  would  always  be  equal  to  unity.  In  our  definition  of 

.T 


B,  it  is  not  true  that  one  of  the  eigen  values  of  BE  would  always  be 
equal  to  unity.  Since  we  know  the  marginal  probability  laws  of  X  and  Y, 


we  need  to  estimate  only  p..'s  by  n../N's.  This  what  motivated  us 

ij  7  ij 


to  define  the  matrix  B  the  way  it  was  defined  above.  However,  in  view  of 


the  above  formula  for  one  could  estimate  x ^  directly  without  having 


to  define  the  matrix  B.  Accordingly,  let 

2 


% 


(nlin22  "  n12n21)  (nlln23  "  n12n21) 

/  '  /  ' 


N  P1P2tJ192 


N  plp2qlq3 


(n12n23  _  n13n22) 


N  PiPoqoq-, 
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We  can  build  a  test  based  on  the  statistic  <2. 

Test  T3:  Reject  the  null  hypothesis  if  <2  1  a* 

We  discuss  the  performance  of  these  three  tests  in  the  case  of  two 
specific  examples  given  below. 

Example  1  p^  =  p2  =  1/2  and  q1  =  q2  =  q3  =  1/3. 

Example  2  Pj  =  1/4,  p2  =  3/4  and  qj  =  1/2,  q2  =  1/4,  q3  =  1/4. 

The  performance  of  T^  and  T3  was  compared  in  detail  in  Subramanyam  and 
Bhaskara  Rao  [12],  We  now  compare  the  performance  of  the  tests  T^,  T2  and  T3 
together  under  the  level  of  significance  a  =  0.01,  0.025,  0.05  and  N  =  15,  20 

A  A  A 

and  25.  The  exact  distribution  of  y,  <  and  k 2  is  evaluated  for  each  of  the 
sample  sizes  N  =  15,  20  and  25  and  the  power  of  the  tests  T^ ,  T2  and  T3  is 

evaluated  at  each  of  the  extreme  point  distributions  of  the  above  examples  using 

these  exact  distributions.  The  graphs*of  these  distributions  are  given  at  the 
end  of  this  section.  Nguyen  and  Sampson  [8]  evaluated  the  powers  of  tests  based 
on  some  other  statistics  at  some  specific  alternative  distributions  by  simulating 
these  distributions. 


*The  authors  wish  to  thank  Ron  Chao  for  his  valuable  help  in  the  computations. 
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Conclusions 

1.  The  power  of  the  test  dominates  the  power  of  the  other  two  tests 
at  the  extreme  point  distribution  P^  in  Example  1  and  Py  in  Example  2 
This  is  not  surprising  as  the  Gamma  Ratio  achieves  the  perfect  value 
unity  under  Pg  in  Example  1  and  Py  in  Example  2. 

2.  On  the  whole  T^  seems  to  perform  well  in  comparison  with  the  other 
two  tests.  Even  under  the  distribution  P^  in  Example  1  and  Py  in 
Example  2,  T^  is  not  overpowered  by  T^. 

3.  Some  extensive  studies  are  needed  to  be  carried  out  to  see  whether 
T^  is  preferrable  to  the  other  two  under  different  sets  of  marginal 


distributions. 
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Portability  Portability  Portability  Portability 


5 .  Inference  on  the  structure  of  dependence 

2 


In  two-way  contingency  tables,  the  x  test  for  independence  has 
been  widely  used.  When  the  test  for  independence  is  rejected,  it  is  of 
interest  to  study  the  structure  of  dependence  between  the  ar+1  rows  and  b+1  columns. 


In  this  section,  we  write  the  matrices  F  =  ( f . .)  and  B  =  (b..)  in  terms 

ij  ij 


of  their  eigenvalues  and  eigenvectors  by  singular  value  decomposition; 

here  fij  *  pij//?tV  bij  =  pi  ■  p11  +  +  pi  ,b+l  ’  Ij  ' 


Pi  .  +  ...+  p  , .  . ,  n .  =  n . .  + 

Klj  a+l,j’  l*  il 


+  n,-  k±i  ancl  n  .=  n.  .  +  ...-  +  n  . 
i  »b+l  •  j  1  j  a+1  ,  j 


Taking  advantage  of  the  above  decomposition,  we  propose  procedures  to 
find  out  as  to  which  of  the  last  t  eigenvalues  of  FF'  are  zero.  The 
distribution  theory  associated  with  the  above  procedures  are  also  in¬ 
vestigated.  Some  aspects  of  the  above  problem  were  discussed  by  O'Neill 
([17],  [18],  [1^]).  The  problem  of  determination  of  the  rank  of 
F  is  discussed  in  a  forthcoming  paper  of  Z.D.  3ai,  P.R.  Krishnaiah  and 
L.C.  Zhao  from  the  point  of  view  of  model  selection  using  an  information 


theoretic  criterion.  The  above  authors  also  established  the  strong  con- 

a+1  b+1 


sistency  of  their  procedure.  In  this  section,  we  use  the  notation  n=  VL  £  n.. 

i=l  j=l  1J 


is  fixed  and  the  marginal  totals  n.  and  n  .  are  random. 

1  O  •  J 


Consider  the  model 


pij  =  Wij 


(5.1) 


i  =  1  ,  2 ,  . . . ,  a  +  1  ,  j  =  1 ,  2 ,  . . . ,  b  +  1  .  Without  losing  general i ty  we 
assume  that  a  <  b.  Under  the  above  model,  it  is  of  interest  to  test  for 


the  structure  of  5^..  From  singular  value  decomposition  of  a  matrix,  it 
is  known  (e.g.,  see  Lancaster  [12])  that 


F  =  £*n*'6  +  T  6  £%*' 

,0.0  o  t,  USUJU 


(5.2) 
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where  6Q  _>  5-j  _>  . ..  _>  6a  are  the  eigenvalues  of  F,  6q  =1,  £*  is  the 

2 

eigenvector  of  FF '  corresponding  to  6^  and  n*  is  the  eigenvector  of 

2 

F ' F  correspondi ng  to  6  .  In  (5.2) 


55 =  (,/hcui *  •••• 

n*  =  ,  .... 

5ol  **  *  “  ?o,a+l 


^a+1  ^u,a+l 
'/qb+lnu,b+l 
=  ]*  nol  = 


)' 

)' 


no  ,b+l 


=  1. 


(5.3) 


We  are  interested  in  finding  out  as  to  how  many  6^'s  are  equal  to  zero. 
This  problem  is  analogous  to  the  problem  of  studying  the  structure  of 
interaction  term  in  two-way  classification  model  with  one  observation  per 
cell.  So,  we  will  briefly  discuss  the  above  model. 

Let 


E(x.  •)  *  u  +  a-  +  B  .  +  A.  .  (5  .4) 

v  ij'  1  J  lj 

for  i  =1,2,  . . . ,  (a  +  1 ) ,  j  =  1 ,  2,  ....  (b  +  1). 


Also  y,  a-,  Sj  and  A . .  respectively  denote  the  general  mean,  i-th  row 

effect,  j-th  column  effect,  and  interaction  in  i-th  row  and  j-th  column. 

Also,  let  A  =  (a..)  and  a  <  b.  We  assume  that  x.  .’s  are  distributed  in- 
*  J  *  J 

2 

dependentally  and  with  variance  a  .  The  problem  of  testing  the  hypothe¬ 
sis  A  =  0  was  first  considered  by  Fisher  and  MacKenzie  [4],  and  later  by 
Williams;  [23],  Tukey  [22]  and  others  when  the  underlying  distribution  is 
normal.  Fisher  and  MacKenzie  [4]  considered  this  problem  using  eigen¬ 


values  of  certain  matrix.  For  a  review  of  the  literature,  the  reader  is 
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re  ferred  to  Krishnaiah  and  Yochmowitz  [10],  Now,  let  EU^)  = 
log  pij. ,  u=  log  c,  ai  =  log  p. ,  8^  =  log  p^.  and  A .  ^  =  log  6.^.  Then 
the  model  (5.4)  can  be  written  as 


Pij  -  CPiqjCij 


(5.5) 


But,  here  we  do  not  assume  that  the  conditions  (5.5)  are  satisfied.  But, 


log  (J>.)  =  log  (£q.)  =  0. 
i  J  J 

We  may  assume  that  c  =  1  ( i  .e . ,  p  =  0) .  We  can  write  (5.6)  as 


(5.6) 


pij  *  c',i<'jexp(nij) 


(5.7) 


and  express  n  =  (n^j)  in  terms  of  its  eigenvalues  and  eigenvectors  using 
spectral  decomposition  of  a  matrix.  Then,  we  can  draw  inference  on  the 
rank  of  n.  This  problem  is  different  from  the  problem  of  drawing  in¬ 
ference  on  the  rank  of  c  except  for  the  special  case  when  the  rank  of 
n  is  zero.  This  special  case  is  equivalent  to  the  statement  that  the 
rank  of  p  is  one.  In  studying  the  interaction  term  in  two-way  classifi¬ 
cation  model,  Tukey  [22]  and  Mandel  [15]  assumed  certain  structures 
on  interaction  term.  We  can  assume  similar  structures  on  the  models 
(5.6)  and  (5.7).  As  far  as  the  models  (5.6)  and  (5.7)  are  concerned, 
they  are  analogous  to  the  well  known  two-way  cl assi fi cati on  model  with 
interaction  and  one  observation  per  cell.  But,  the  problems  of  estima¬ 
tion  and  distributions  of  test  statistics  are  of  different  nature.  In 
general  we  may  also  consider  models  of  the  form 

Pij =  f(“c  “r  "lj> 


« .  i  ^  I  ^  N.  "V  V  •  .  * 
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where  f ( • )  is  a  suitably  chosen  function  of  a-j  ,  3  ^ ,  and  n^..  For 
example  f(x,y,z)  =  f-j ( x) f2(y) f3(z) .  As  another  possibility,  we  may 
also  consider  the  model 


U  +  a.  +  B  .  +  y .  . 
l  J  l  J 


(5,8) 


as  in  two-way  classification  with  interaction. 

Goodman  [61  discussed  the  model  (5.7)  when  n  is  written  in 
terms  of  its  eigenvalues  and  ei genvectors .  O'Neil  ([171,  [10],  T 1 9] ) 
discussed  some  aspects  of  the  asymptotic  distribution  theory  associated  with 
finding  the  rank  of  the  matrix  c .  In  this  section,  we  propose  various 
test  procedures  for  determination  of  the  rank  of  t  and  investigate  some 
problems  on  the  asymptotic  distributions  of  the  test  statistics. 

We  now  discuss  the  problem  of  testing  for  the  rank  of  the  matrix  t . 

If  we  know  in  advance  the  rank  of  t,  we  can  use  that  knowledge  in  esti¬ 
mating  the  unknown  parameters  more  accurately.  For  example,  the  maximum 

likelihood  estimates  of  p..'s  when  the  rank  of  F  is  one  are  not  the  same 

J 

as  when  the  rank  of  F  is  greater  than  one. 

Now,  let  B  =  (b^  .)  where 


n.  .//nT 
i  j  t 


(5.9) 


Then,  using  the  spectral  decomposition  of  B,  we  have 


A  A  A 


B  =  +  T  6  4*0*' 

~o~°  0  u=i  u~u~u 


(5.10) 


where  5=1, 
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"u  '  ^n.i/n)  ^  ’  ‘  *  *  ,(n*b+l,/n^ 


1/2 


nu,b+l) ' 


?ol 


*  ‘  =  ?o,a+l  ”  1  '  nol 


o  ,b+l 


=  U  60  =  1 


Also,  6q>^6-j  >_  ...  ^.<Sa  are  the  eigenvalues  of  B,  £*i  is  the  eigen- 


~2 

vector  of  BB  corresponding  to  6u  and  n*  is  the  eigenvector  of  B'B 


„  2  ?  a 

corresponding  to  6^.  Now,  let  H.:6(r  =  0  (i  =1,2,  ...,  a)  and  H  =  /J.H.. 

We  can  use  ^(o,*  ....  <S  )  to  test  H  where  ^(*)  is  a  suitable  function  of 

^  >  a 

'2  "2  "2  “2  ~ 2 

5,,  ....  6  .  For  example,  we  can  use  6,  +  ...  +  6  or  5,  as  test  statistics. 

la  I  d  1 

*2  "2  2 
Here  we  note  that  n(6^  +  ...  +  6&)  is  equivalent  to  xQ  where 


X2Q  «  -  (ni.n.j/n)>2/ni.n.j-  (5J2> 

1 J 

When  H-j  is  true,  O'Neil  [17]  showed  that  the  joint  distribution  of 
“2  "2 

no.,  ...,  n6  is  the  same  as  the  joint  distribution  of  the  eigenvalues 
of  the  central  Wishart  matrix  with  a b  degrees  of  freedom.  Percentage 
points  of  the  largest  eigenvalue  of  the  central  Wishart  matrix  are  given 
in  Krishnaiah  [8], 

We  have  discussed  before  some  procedures  to  test  for  the  overall 
2  2 

hypothesis  6,  =  ...  =  6a  =  0.  We  will  now  discuss  procedures  for  testing 
the  subhypotheses  H  when  H-|  is  rejected.  The  hypothesis  is  the  same 

as  the  hypothesis  that  the  rank  of  F  is  t.  We  will  first  consider  the 

'2  ~2 

test  procedure  based  upon  T,  where  T  =  6  +  ..  +  5  .  In  this  pro- 

I  q  q  a 

cedure,  we  accept  or  reject  according  as 


where 


T1  >  Cla 


(5.13) 


Pr-T!  1  clalH!‘j  =  O-*)- 


(5.14) 


If  is  rejected,  the  hypothesis  Hq  is  accepted  or  rejected  according 


T  <  c,  . 

q  >  la 


(5.15) 


Starting  with  the  test  based  upon  we  can  also  draw  inference  on  test¬ 
ing  the  hypothesis  H^.  as  described  below;  here  .  denotes  the  hypothesis 

2 

p.  .  =  p.q.  for  given  values  i,j.  The  x  test  statistic  can  be  written  as 


where 


x  =  z'z 


(5.16) 


-  *  ~  (  7  7 

\  ^-1  1  >  •  o  •  » 


11  ’*°*,21  ,b+l  ’ ' *  * ’za+l  ,1  *  *  *  * ,Za+l ,b+l 


But  xQ  3  rnax(c'z)"  when  the  maximum  is  taken  over  all  non-null  c  subject 
to  the  restriction  that  c'c  =  1,  So,  when  H-|  is  rejected,  we  can  test 
the  subhypotheses  H..  as  follows.  We  accept  or  reject  H..  against  two- 
sided  alternatives  according  as 


z .  .  <  c,  . 

1 .1  >  «  a 


(5.17) 


t  u 

We  can  test  the  hypothesis  •'  T  H.  .  as  follows.  We  accept  or  reject 

i =1  j=l  1J 

the  above  hypothesis  against  two-sided  alternatives  according  as 

t  u  - 

v  v  ,2  < 

i-1  j  =  l  1j  >  CH 


t  u 

The  hypothesis  0  OH.,  implies  the  hypothesis 

j  =  l  j  =  l  1J 
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We  will  now  discuss  the  problem  of  testing  the  hypotheses  H..  .  against 

the  alternatives  A.,  simultaneously  where  A..:p..  >  p.q..  We  accept 
i  j  i  j  i  J 

or  reject  H.  .  against  A. ^  according  as 

z.  .  <  c0  (5.18) 

ij  >  2a 

where 

P[max  z..  <_  C2alH-|]=  (1-a)  (5.19) 

and  max  2. .  denotes  the  maximum  of  the  elements  of  2,  The  asymptotic  joint 
distribution  of  the  lelments  of  2  is  a  singular  multivariate  normal  distribution. 
But,  bounds  on  the  critical  values  c2a  can  be  obtained  by  using  Poincare's 
formula.  Sirrrilarily,  we  can  propose  procedures  to  test  hypotheses  H.  . 

J 

against  At^  where  ATj  :pi  .  <  p.q^. 

*2 

We  now  discuss  the  test  based  upon  p-j .  We  accep.t  or  reject 
according  as 

;2<c2a  (B.20) 

where 

p[p2  >  c2alH1J  =  (1-a).  (5.21) 

If  H-|  is  rejected,  we  accept  or  reject  H.  according  as 

p2.  >  c2  .  (5.22) 

We  will  now  derive  the  asymptotic  nonnull  distributions  of  certain 
functions  of  no.,..., no  .  The  following  lemmas  are  needed  in  the  sequel. 

Lemma  5 . 1  Let  U :  p  x  p  be  a  symmetric  matrix  which  can  be  expressed  as 

U  =  A  +  e  l/1  ^  +  2 


+  . . . 


(5.23) 


when  A  =  diag( , . . . ,A  ) .  Also,  let 


Qi  ■*"•  •  •  to  i 


.  =  l  ^  A  =  e 

!-,+.» .+Qa  a 


for  a  =  l,2,...,r,  qQ  =  0  and  qr  =  a.  In  addition,  let  > 
denote  the  eigenvalues  of  U.  Then 


i  =  e 

a.  a 


+  £i">  +e¥2>  ♦ ... 


(5.24) 


*  1  "p 


(5.25 


where 


!<’)  =  — tr  U<>> 
q  act 
a 


J2)  ^  —Ur  |u^  +  l  Q~] 

Q  .  rsrv  r^K  rv  K  ft  rv 


aa  a6  aB  Ba 

Bra 


0  0  =  0  -  0Q 
aB  a  B 


r  ij(1)  y(i)  u(d 
un  u12  •••  ulr 


U  = 


y(i)  y(i ) 
U21  U22 


y(i)  y(i) 

Url  Ur2 


..  U 


(i) 

2r 


3(i) 

rr 


The  above  lemma  is  implicit  in  Kato  [24],  It  is  also  proved  in 
Fujikoshi  [5]  by  following  the  same  lines  as  in  Lawley  [13], 


Lemrna  5.2  Let  U,  , . . .  ,t) ,  ( i  =1 ,2 ,. . .  ,k) ,  be  analytic  functions  of 

- . -  1  I  a 

Z,,...,?.  around  A,,...,A  and  let  Ms  have  multiplicities  as  in  (5.24) 

13  13  1 


We  assume  that 


for  e  J  ^  *  J*2  ^  ^  g  *  J  3  ^  ^  X  ”  (\j»#««>x^)  3 nd  J ^ 

denotes  the  set  of  integers  q^  +...+  q^  +  l,...,q^  +...+  qa.  Then 

Li  "  ^  <VV-“.2a)  ~  VAl . Aa>> 


/S  a-1  iQ  « 


aB  a0  w0a 


+  I  n»la8tru“,trU«>+...J  (5.27) 

a  0 


for  i  =  l,?,...,k.  Let  H  and  K  be  orthooonal  matrices  and  let 


R1  =  H'BK. 


If  we  choose  the  first  columns  of  H  and  K  such  that 
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0  R 


Similarly,  let  =  H*'FK*  where  the  first  columns  of  H*  and  K*  are  given  by 


hiO  =  pi 


*  i 

kj0  =  qj 


./l  o 


0  fi 


/I  0  \ 


1  0 


RtRJ  =:'  ■  ,  fi,fi]  = 

1  1  o  RR'  1  ,0  nn' 


A  2  A  2 

It  is  known  (see  O'Neill  [17]  )  that  p1  >_...>  og  are  the  eigenvalues  of  RR' 

o  2 

whereas  oT  >...>  are  the  eigenvalues  of  fin'.  Next,  let 

i.  —  —  a 


X  =  ft  (R-fi)  =  (xiJ) 


(5.28) 


where  R  =  (r..)  and  fi  =  ( co . . )  -  Then  x,.'s  are  known  (see  O'Neill  [17])  to 
1 J  1 J  *  J 

be  asymptotically  distributed  as  multivariate  normal  with  mean  vector  0  and 
the  elements  of  the  covariance  matrix  are  given  by 


Cov(xij.  ,xkJ ,)  -  ai j . kg,  (say). 


(5.29) 


Using  (5.28) ,  we  obtain 


Now,  let 

Then, 

where  V  : 
=  M'XX'M 
and  's 


But, 


RR'  =  nn'  +  n"*(fiX‘  +  XQ')  +  n_1XX' 


(5.30) 


M:  axa  and  L:  bxb  be  such  that 


I'nL  =  (diag(p1,...,Pa)|0)  =  Dp 


+  jl  vw  ;iv^v 


S  =  M'RR'M  =  V  +  —  +  ~r  V' 

/n  n 


(5.31) 


D  d:,  V(1)  =  (V^)  =  dj1  +  ZD*.  Z  =  (z..)  =  M'XL ,  and  V(2)  =  (/2) 

*2 


p  p  '  '  a&  '  ~p~  P'  '  ij'  *  “6 

(1)  M)  -  .  -  —-‘2.-2 


,  Here  V  „  and  V'„  are  of  order  q  xq  Now,  let  i.  -  p.,  X.  =  p. 

aB  a8  a  8  i  i  i  i 

have  multiplicities  as  in  (5.24).  Then,  using  (5.26),  we  obtain 


L.j  =  </n  {T.  ( n..  ,oo  •  ,£p)  -  T^.  ( X j ,. . .  ,Xp) } 


=  V  a-  tr  V 

L  -  la  aa 


a=l 


(1)  1  f  T  a  tr{V^  +  y 

+  -L.  L  i  aiatri  aa  A  oB  ap  V 
aa  /—  a-i  Bra 


/n 


+  R  tr  V(1)  tr  vj^] 

LL  laB  aa  ftR 


86 


+  terms  of  higher  order. 


(5.32) 


,(D 


2plzll 


(p2Z12  +  P1Z21)  ’*•  (paZla  +  plZa! 


2p2z22 


(plz21  +  P2Z12) 

(PiZ  ,  +  P„ziJ  (p2Za2  ■  a 


(paZ2a  +  p2Za2} 


P-^0_)  ...  ^p 


'2a' 


a  aa 


V 


(2) 


M^XX'M1  M^XX'M2  ...  M^XX'Mr 
M^XX'M2  ...  M^XX'Mr 


M'XX'M.  M’XX'M,  ...  M'XX'M 
^  a  r  *»  r  r 


where  M  =  (M^,..,,Mr)  and  NK  is  of  order  axq...  So, 


r 

'i  "  2a=iaia9a(Zqi+---+Vi+1»qi+--+Vi+1  +“‘+  zq1+...+qa’qi+---+qa) 


=  5i?0 


(5.35) 


where 


(2ailel^ 


,...,2a1rery 


(5.36) 


and  z'0  =  (z^,...  ,zaa).  The  asymptotic  distribution  of  B'z  is  multivariate 

normal  with  mean  vector  0  and  covariance  matrix  B'z  B  where  t  is  the 

•» 

covariance  matrix  of  z  and  B  =  (bj,...,b^).  We  can  summarize  the  above  results 
as  follows: 

Theorem  5.1  We  assume  that  p.'s  have  multiplicities  as  given  below: 


ql+. . *+qa_i 


.=  p 


q,+. ..+q 

1  a 


(5.37) 


for*  a  =  1,2,... ,r  where  qQ  =  0,  q^  +...+qr  =  a.  Also,  let  L1,...,L)<  be  functions  of 
~2  -2 

°l**,*»pa  satl" sfy' n9  the  assumptions  (5.25).  Then,  as  n  ->  ®,  the  joint  distri¬ 
bution  of  L1,...,Lk  is  multivariate  normal  with  mean  vector  0  and  covariance 
★  * 

matrix  B  z  B  where  i  is  the  covariance  matrix  of  z,  B  =  (b1f,#>1b.)  and  b.'s 


are  defined  by  (5.36). 


ft 


Now,  let  V  a.  trV'1^  =  0  for  each  i.  Then, 

lei  aa 

a=l 

-  n{T..  (s,j  ,. . .  ,£p)  -  ( Xj ,Xp) } 

/Vi' 

-  i(vx . vr)  A/  | 1 

V' 

where 

Ai  =  ^laB* 

vi  =  2^pizn  +,*,+  Pq1zq1q1^ 

V2  =  2{pq1+lZq1+l,q1+l  +  ***  +  ^^q^q^ 

v  =  2  ( n  7  + 

r  ^VPq1+..o+qr_1+rq1+...+qr_1+l,q1+...+qr_1+l  *** 

q^+.  .<.+qr  q  j"*-*  •  •  "^qp »q  j-*"*  •  •  "^q j,)  < 

Since  p/s  have  multiplicities,  we  can  write  v  =  ,v  )’ as  v  =  Ez0  where 

E  =  (5’-**’er)j 

e.l  •  29i<iij- 0 . 0) 

eX  =  2e?(0...  ,0,1’  ,0,...,0) 

•mC.  L  2 

?r  ■  29r'° . °’^) 

As  n  •*■»>,  v  is  distributed  as  a  multivariate  normal  with  mean  vector  0  and  co- 

if  if  if 

variance  matrix  Ez  E'.  So,  the  joint  distribution  of  is  the  same 

as  that  of  correlated  quadratic  forms  discussed  in  Khatri ,  Krishnaiah  and  Sen  [7] 
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